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Abstract 

The scattering potential of the defects on Ge(001) surfaces is investigated by first-principles 
methods. The standing wave in the spatial map of the local density of states obtained by wave 
function matching is compared to the image of the differential conductance measured by scanning 
tunneling spectroscopy. The period of the standing wave and its phase shift agree with those in the 
experiment. It is found that the scattering potential becomes a barrier when the electronegativity 
of the upper atom of the dimer is larger than that of the lower atom, while it acts as a well in the 
opposite case. 

PACS numbers: 73.20.Hb, 72.10.Fk, 82.20.Wt 



7 With the progress of new techniques for the atomic-scale manipulation and modification 

8 of materials, there is considerable interest in the electron scattering properties of nanos- 

9 tructures both experimentally and theoretically. Scanning tunneling microscopy (STM) |l| 

10 and mechanically controllable break junctions [2J have been extensively used to study the 

11 conductance of atomic-scale systems suspended between two electrodes so far, and today it 

12 is already well known that the conductance of a single strand of atomic wire is quantized in 

13 the unit of 2e 2 /h, where e is the electron charge and h is Planck's constant {3)]. Although 

14 these techniques are powerful tools to understand electron transport phenomena of such 

15 minute systems, it has not been easy to discuss the contribution of local chemical bonds 

16 to the electron scattering because these techniques only measure the current between two 

17 electrodes with an applied bias voltage. On the other hand, the spatial maps of the local 
is density of states (LDOS) obtained by scanning tunneling spectroscopy (STS) can provide 

19 the images of standing waves, which give important information about the dispersion rela- 

20 tion of surface states [J, |5|] as well as about the electron scattering process at the potential 

21 barrier 4j. 

22 Tomatsu et al. demonstrated the standing wave around oppositely buckled Ge dimers 

23 on a Ge(001) surface and discussed the difference in scattering potential between the Si- 

n q 

24 Ge and Sn-Ge dimers [y, [7j. They also carried out first-principles calculations based on 

25 the density functional theory (DFT) p] using periodic supercells to identify the surface 

26 atomic and electronic structures observed in STM images and STS spectra. However, the 

27 standing waves emerge as a result of the difference in the coefficients between incident and 

28 reflected waves of scattering wave functions, which cannot be calculated by models under 

29 periodic boundary conditions along the scattering direction. Thus, it is not straightforward 

30 to determine the scattering potential using the periodic models. 

31 The scattering properties of atomic-scale structures are studied by DFT calculations 

32 combined with nonequilibrium Green's functions |9[ or wave function matching methods 



33 Il2l|. In these approaches, the scatterers are sandwiched between electrodes that extended 

34 semi-infinitely and are connected to bulks; thus, wave functions extending over the entire 

35 system can be correctly described as a scattering state distributed as a result of the existence 

36 of scatterers when an electron comes from infinitely deep inside the electrode. Although these 

37 approaches have been used to investigate the transport properties of atoms, molecules, and 

38 thin films suspended by two electrodes so as to model the experiments using the STM and 



39 mechanically controllable break junctions, to the best of my knowledge, no studies have 

40 attempted to obtain an understanding of the relationship between local chemical bonds and 

41 their scattering properties by examining the scattering potential in first-principles transport 

42 calculations. 

43 In this Letter, I apply the first-principles transport calculation to the investigation of the 

44 standing waves due to scattering at Si-Ge and Sn-Ge dimers on a Ge(001) surface. The real- 

45 space procedure of the DFT calculation combined with the overbridging-boundary matching 

46 (OBM) method provides the reflection and transmission coefficients of the incident electrons, 

47 which enable us to determine the scattering potential and the phase shift of the standing 

48 waves. The calculated phase shift of the standing waves is in good agreement with that 

49 measured in experiments Q. The defects act as a barrier when a Si (Sn) atom exists at the 
so lower (upper) position of the dimer, while they behave as a well in the case of the other 

51 way around. It is found that the scattering potential is related to the stabilization of the ir 

52 bands of the Ge(001) surface due to the Jahn- Teller effect of the dimers and the difference 

53 in electronegativity between Ge and the impurity atom. 

54 The first-principles calculation method used here to obtain the electron scattering prop- 



55 erties is based on the real-space finite-difference approach 



12 



-Il4j. which enables us to de- 



56 termine the self-consistent electronic ground state with a high degree of accuracy by a 



57 timesaving double-grid technique 



12 



14j . Moreover, the real-space calculations eliminate 
53 the serious drawbacks of the conventional plane-wave approach, such as its inability to de- 

59 scribe nonperiodic systems accurately. The scattering properties of Ge-Si and Ge-Sn dimers 

60 on a Ge(001) surface suspended between Ge(001)-(2x2) surfaces are examined by making 
ei use of this advantage of the combination of the real-space finite-difference approach and the 
62 OBM method Lli LA 



HI 

63 Figure [1] shows the computational model, where the scattering region is connected to 

64 the Ge(001)-(2x2) surfaces. The scattering region consists of six Ge atom layers with a 

65 (2x6) lateral cell, and the topmost atoms are buckled so as to form a surface reconstructed 

66 structure, while the bottommost layers are terminated by hydrogen atoms. The dimensions 

67 of the cells are L x = 23.92 A, L y = 7.97 A, and L z = 17.94 A for the scattering region, 
es and L x = 7.97 A, L y = 7.97 A, and L z = 17.94 A for the Ge(001)-(2x2) surface electrode 

69 region, where L x and L y are the lengths in the x- and y-directions parallel to the surface, 

70 respectively, and L z is the length in the ^-direction. The Ge-Ge dimer at the center of the 



71 scattering region is replaced by Ge-Si or Ge-Sn dimer as a defect. When the Si and Ge atoms 

72 are located at the lower and upper sides of the dimer, respectively, the dimer is referred as 

73 an SiL dimer. Other dimers are named in a similar manner. During the optimization of 

74 the atomic configurations in both the electrode and scattering regions, all atoms except 

75 for those in the bottommost layers are relaxed until all the force components drop below 

76 0.05 eV/A by imposing the periodic boundary condition in the x-, y-, and ^-directions. 

77 The grid spacing is set at ~ 0.20 A and A;-space integrations are performed with 4 x 4fc 



78 points in the irreducible wedge of a two-dimensional Brillouin zone of the Ge(00T 



•(2x2) 



79 surface. Exchange correlation effects are treated by a local density approximation [15| of the 
DFT, and the projector augmented wave method [3] is used to describe the electron-ion 



80 

8i interaction 



82 In the OBM calculation per 



ormed to obtain the scattering wave functions, the norm 



83 conserving pseudopotentials \Vj§ of Troullier and Martins [18] are employed instead of the 

84 projector augmented wave method. To determine the Kohn-Sham effective potential, a 
as supercell is used under a periodic boundary condition, and then the scattering wave functions 
se are computed under the semi-infinite boundary condition obtained non-self-consistently. It 
87 has been reported that this procedure is just as accurate in the linear response regime but 
ss significantly more efficient than performing computations self-consistently on a scattering- 

89 wave basis 

90 Figure [2] shows the charge density distribution of the scattering wave functions for the 

91 electrons propagating from both the left and right electrodes with an energy of -Ej? + 0.55 eV, 

92 which corresponds to the spatial image of the differential conductance in the STS spectrum 

93 0]. Here, Ep is the Fermi energy. The standing waves on the xy plane, which is 2.7 A 

94 from the upper Ge atom of the Ge(001) surface, are plotted. It is found that the standing 

95 wave emerges in all the models. Since the 7r* band of the Ge(001)-(2x2) surface crosses the 

96 energy of Ep + 0.55 eV at k x = 0.454 x n/L x , the period of the standing waves is 2k x (=17.6 

97 A), which is slightly shorter than that observed in the experiment (~ 22 A). Taking into 

98 account the fact that the local density approximation underestimates the energy band gap 

99 and that the period of the standing wave decreases as the sample bias increases in the spatial 

100 image of the STS spectrum in Ref. Q, these results agree with those of the experiments. 

101 The line profiles of the standing waves along the dimer row including the impurity atom 

102 are depicted in Fig. [31 To demonstrate the period and phase shift of the standing waves 



103 clearly, the amplitude of the waves is fitted by 



a(x) = A cos(2k x + (/)), (1) 

104 where A and (f> are the amplitude and phase shift of the standing waves, respectively. The 

105 fitting is carried out using the amplitude of the standing waves on the upper atom of the 
loe dimer (indicated by "•" in Fig. [2]). Table [I] shows the amplitude and phase shift of the 
107 standing waves. The amplitude of the SnL (SiU) dimer is larger (smaller) than that of the 
los SiL dimer, and the phase shifts of the SiU and SnL dimers are negative, while that of the 
109 SiL dimer is negative. These findings consistent with experimental results 

no Let us discuss the scattering property of the defects in more detail. The OBM method 
in provides the reflection coefficients of the scattering waves, which give us information on the 
n2 scattering potential upon using a one- dimensional free-electron-like model. Table [III shows 
n3 the calculated reflection coefficients and the reflection probabilities. The barrier height V 
H4 and length a of the scattering potential are determined by 

ref = (k 2 - K 2 )(l - e 4iKa )e 2ika 

° (k + K) 2 - (k-Kye 4 ***' [ ' 

H5 which is derived from the penetration of electrons into a one- dimensional square potential 
H6 barrier on the basis of quantum mechanics. Here, k = v/h and K = y (v 2 — 2mV)/h, where 
ii7 h is the reduced Planck's constant, m is the electron mass, and v is the group velocity of 
us the incident wave. V and a are fit so that the reflection coefficients correspond to those 
n9 obtained by the OBM method. The calculated V and a are shown in Table II. The SiU 

120 and SnL dimers behave as a potential barrier, whereas the SiL and SnU dimers behave as a 

121 potential well. The signs of the scattering potential height except for that of the SnU dimer 

n 

122 agree with those reported by Tomatsu et al. |7| , where the potentials are estimated by fitting 

123 the line profile of the STS spectrum and/or by DFT calculation using periodic supercells. 

124 However, the potential of the SnU dimer, which is not observed in the experiment, is not in 

125 agreement. 

126 It is known that an electron conducts through the ir* band of the Ge-Ge dimer on the 



127 Ge(001) surface. According to the explanation by Haneman 20j, the sp 3 bonding state 

128 is stabilized as the bond angle around the atom of the dimer decreases. In the case of a 

129 Ge(001) surface, the n band accumulating around the upper atom of the dimer is occupied 

130 after the surface reconstruction, as shown in Fig. H^a). When the electronegativity of the 



131 upper atom is larger than that of the lower atom, the energy of the it (jt*) band further 

132 decreases (increases) owing to the low (high) electronegativity of the lower (upper) atom 

133 [Figs. @Jb) and@Jc)]. Since the Mulliken electronegativities of Si, Ge, and Sn are 2.0, 1.9, 

134 and 1.8, respectively, the energy of the 7r* bands of the SiU and SnL dimers is higher than 

135 that of the SiL and SnU dimers [Figs. @Jb) and@|c)]. Thus, the scattering potential of the 
we SiU and SnL dimers acts as a barrier for the conducting electrons, whereas that of the SiL 
13? and SnU dimers acts as a well. In addition, because the electrons conduct through the states 
us accumulating at the lower atom of the dimer, the length of the scattering potential barrier 
139 increases long when the impurity atom exists at the lower side of the dimer. In Ref. Q, the 
wo projected charge density distribution of the tt band of the SiL dimer is high at the defect, 
wi while it is low in the case of the SnL dimer. Moreover, the energy gap between the tt and 

142 ?t* bands, e n * — e n , of the SnL dimer is larger than that of the SiL dimer, which supports 

143 the explanation. 

144 In summary, the scattering potential of the defects on Ge(001) surface is investigated by 

145 first-principles calculation. The phase shifts of the standing waves due to the defects are in 
U6 agreement with those obtained by experiments |7j ■ By calculating the reflection coefficients 
147 of the scattering wave functions, it was found that the scattering potentials of the SiL and 
us SnU dimers act as a well, while those of the SiU and SnL dimers behave as a barrier. 

149 This characteristic is interpreted in terms of the electronegativity of the defects; when the 

150 electronegativity of the upper site of the dimer is large, the energy gap between the n and 

151 7r* bands increases, resulting in the generation of the potential barrier for the it* electrons. 

152 This explanation is also consistent with the LDOS in the vicinity of the dimers and with the 

153 value of En* — e n computed using a conventional periodic supercell. In addition, the results 

154 demonstrate the applicability of first-principles transport calculations to the investigation 

155 of the scattering potential and phase shift due to defects on surfaces in collaboration with 
we STM and STS measurements. 
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FIG. 1. Schematic image of computational model. Dark (blue) spheres are Ge atoms and light 
(yellow) spheres represent the dimer replaced by a Ge-Si or Ge-Sn dimer. 
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FIG. 2. Surface atomic structure (a) and spatial images of standing waves for (b) SiL, (c) SiU, (d) 
SnL, and (e) SnU dimers. Each contour represents a density of 0.496 x 10~ 5 e/A 3 /eV/spin higher 
or lower than that of the adjacent contours. 
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FIG. 3. Line profiles along thick solid lines in Fig. [2] for (a) SiL, (b) SiU, (c) SnL, and (d) SnU 
dimers. Solid and dashed curves represent the standing waves and the curves fitted by Eq. ([1]) 
using the data indicated by bullet points, respectively. 
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207 TABLES 



208 



TABLE I. Amplitude and phase shift of the standing waves around the Ge-Si and Ge-Sn dimers. 
Model Amplitude A (e/A 3 /eV/spin) Phase shift (rad) 

SiL 3.908xl(T 7 0.2217T 

SiU 1.692xl(r 7 -0.602tt 

SnL 7.169xl(T 7 -0.650tt 

SnU 1.031 xl(T 7 0.1427T 



16 



TABLE II. Reflection coefficients obtained by OBM method, barrier heights, and barrier lengths 
calculated from the reflection coefficients using Eq. ([2]). 
Model Coefficient Height V (V) Length a (A) 

SiL 0.0519 + 0.0565i -0.801 5.203 

SiU -0.0044 - 0.0270i 4.444 3.518 

SnL -0.0293 - 0.1092i 1.394 4.858 

SnU 0.0138 + 0.0097i -0.152 4.378 
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